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Abstract 

The Kekule structure of the ground state of (5,5) armchair carbon nanotube is revealed by semiempirical molecular orbital calculations. This 
structure has bonds with two different bond lengths, differing by 0.003 nm. The ground state has tripled (compared to undistorted case) 
translational period due to Peierls distortions. Two first order structural phase transitions controlled by the tension are predicted at zero 
temperature. These transitions correspond to 5% and 13% elongations of a uniaxially deformed (5,5) nanotube. The narrow gap semiconductor 
to metal transition is predicted at 5% elongation of the nanotube. 



1. Introduction 

The studies of electronic and elastic properties of carbon 
nanotubes (CNTs) are actual in connection with perspectives 
of their applications in nanoelectronic devices [1] and in com- 
posite materials [2], and are also of fundamental interest, par- 
ticularly for physics of phase transitions. For example, the su- 
perconductivity [3], the commensurate-incommensurate phase 
transition in double-walled nanotubes [4] and the spontaneous 
symmetry breaking with formation of corrugations along nano- 
tube axis in ultrasmall CNTs (0.4 nm in diameter) [5] have 
been considered. 

Within a framework of the tight-binding model the electronic 
properties of CNTs are determined by their diameter and chi- 
rality [6] and armchair (n, n) CNTs are metallic. The first clear 
experimental evidence of the chirality dependence of electronic 
properties of CNTs have been gained in STM/STS measure- 
ments [7,8]. However such kind of measurements can not re- 
solve any bond length difference in CNTs and now [1] elec- 
tronic properties are conventionally considered for CNTs with 
all equal graphitic C-C bond lengths. 

The possibility of Peierls transition in carbon nanotubes was 
first considered in [9]. As a result of this transition, armchair 
CNTs become semiconducting at low temperature and Peierls 
distortions lead to Kekule structure (see Fig. 1(a)) with triple 
translational period Lf (three times more hexagons in trans- 
lational unit cell). In early works the Peierls gap and Peierls 
transition temperature were estimated in the framework of Su- 
Schrieffer-Heeger model [10, 11] and by means of phenomeno- 
logical electron-phonon coupling parameters [12,13]. For a re- 
cent review on this problem see [14]. 
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By now the Kekule structure of the ground state was found 
only for finite length nanotubes [15-17]. Density functional 
theory (DFT) calculations for finite length (5,5) CNT [15,16] 
give 60 atom periodicity of physical properties on the length 
of CNT segment which is consistent with Kekule structure for 
infinite armchair CNTs. Moreover, X-ray crystallographic anal- 
ysis of recently chemically synthesized short (5,5) CNTs [17] 
show the Kekule bond length alternation pattern for their geo- 
metrical structure, which is in good agreement with DFT and 
PM3 calculations also performed in [17]. 

In the present Letter we consider the possibility of struc- 
tural phase transitions connected with spontaneous symmetry 
breaking and controlled by uniaxial deformation of CNTs on 
the example of the infinite armchair (5,5) CNT. We present the 
PM3 [18] calculations of the internal energy and geometrical 
structure of a uniaxially deformed infinite armchair (5,5) CNT 
up to 17% elongation. For calculated geometrical structure we 
use simple tight-binding (Hiickel) [19] calculations (taking into 
account only the first nearest-neighbour interactions) to gain 
qualitative insight in the electronic energy band structure of a 
uniaxially deformed armchair (5,5) CNT. 

2. Methodology 

The MOPAC-2006 code [20] was used in PM3 calculations of 
the internal energy and the structure of the (5,5) CNT. The ade- 
quacy of the PM3 parametrization of the Hamiltonian has been 
already demonstrated [21] by the calculation of bond lengths 
of the Ceo fullerene with Ih symmetry: the calculated values 
of the bond lengths coincide with the measured ones [22] at 
the level of experimental accuracy of 10"'' nm at liquid helium 
temperature. The PM3 parameterization gives also correct C-C 
bond length in graphite [23]. Comparison of different types of 
semiempirical calculations [24] shows that the PM3 parameter- 
ization gives the best accuracy for structure and internal energy 
of carbon nanostructures. 
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Fig. 1 . Atomic structures of (5,5) carbon nanotube corresponding to: (a) three energetically equivalent structures of phase A (the ground state of the undeformed 
nanotube); (b) phase B of deformed nanotube; (c) quinoid structure of the phase C of deformed nanotube and also of the transition state of the phase A. 
Dotted lines in (a), (b) and (c) are the primitive unit cells of the carbon nanotube for the phases A, B and C, respectively. The coordinate system x, y is 
used for the calculation of the Hamiltonian matrix (2). Translational periods = /2 = = ao\/3 = Lt, if all bond lengths are equal to ao. 



Since the Kekule structure of the infinite armchair CNTs has 
the triple translational period, the proper choice of the length 
of a computational cell along the CNT axis is necessary for 
adequate calculations of their geometrical structure. Namely, 
if one chooses the number of atoms in the computational cell 
of infinite (n, n) CNT which is not divisible on 12n, then the 
Kekule structure is excluded from consideration. This is the 
case in many quantum-chemical calculations of infinite arm- 
chair CNTs (see, e.g., [23,25,26]). For the best of our knowl- 
edge (see also [15]) there is only one quantum-chemical calcu- 
lation [27] of infinite armchair CNT with properly chosen com- 
putational cell for considering the Kekule structure. This early 
version of Complete Neglect of Differential Overlap semiem- 
pirical CNDO/2 calculation gives the quinoid type structure (see 
Fig. 1(c)) for the ground state of the (5,5) CNT being slightly 
lower in energy than the Kekule structure. On the other hand, 
our [28] semiempirical NDDO-type (Neglect of Diatomic Dif- 
ferential Overlap) calculations with the PM3 [18] parameteri- 
zation (as discussed above it is very well suited for the calcu- 
lations of pure carbon systems) give the Kekule structure for 
the ground state of the (5,5) CNT. 

In calculations of the ground state geometrical structure of 
the infinite undeformed (5,5) CNT we used computational cells 
with 60, 120 and 240 atoms in order to check the convergence 
versus the size of the computational cell. Born- von Karman 
periodic boundary conditions along the nanotube axis and the 



full geometry optimization (without any symmetry constraints) 
including the computational cell length were used. 

3. Ground state structure of undeformed (5,5) CNT 

The calculated Kekule structure of the undeformed (5,5) 
CNT ground state at zero temperature is shown in Fig. 1(a). 
Geometrically, this structure is possible for all armchair (n, n) 
CNTs. Three energetically equivalent atomic structures corre- 
sponding to the ground state are shown in Fig. 1(a). The tran- 
sition state between these minima is shown in Fig. 1(c). The 
calculated energy difference between the ground (with tripled 
period) and transition states (with the initial period) of the (5,5) 
CNT is 3 meV per a carbon atom. The bond lengths have the 
values a = 0.1405, b = 0.1433, c = 0.1406, d = 0.1434 nm 
for the ground state and a = 0.1420, b — 0.1426 nm for the 
transition state. Note that the bond length difference between 
short (a and c) and long bonds (b and d) of the ground state is 6 
times greater than the difference between the bonds of the tran- 
sition state. As opposed to the transition state structure with all 
atoms lying on cylindrical surface, the Kekule structure repre- 
sents the alternating rings of atoms with different radii, i?i — 
0.341 nm for rings with pi and qi atoms and i?2 = 0.344 nm 
for rings with p2 and q2 and with and ^3 atoms. Note that 
the same pattern (Kekule structure) with bond lengths close 
(with relative error < 1 %) to our calculations is revealed by 




DFT method for finite length hydrogen terminated (5,5) CNT 
C160H20 [15]. 

4. Deformational structural phase transitions 

Since the translational period of deformed nanotube a priori 
is not known, we use computational cell of 120 carbon atoms 
to match both possibilities of tripling and doubling of the trans- 
lational period of deformed CNT. The structure and internal 
energy of the (5,5) nanotube are calculated for different elon- 
gations (all other geometrical parameters have been optimised) 
£ = {I — leq)/leq, whcrc I is the length of the computational 
cell, and Zeq is the length of this cell, corresponding to the elon- 
gation e = 0. The dependences of the internal energy U and 
the non-equivalent bond lengths with respect to e are shown in 
Fig. 2(a) and Fig. 2(b), respectively. 

The following phases with different symmetries of the struc- 
ture have been found for the uniaxially deformed (5,5) CNT: 
the phase A with the Kekule structure and the triple transla- 
tional period Lf at small elongations e, the phase B with the 
double translational period Lf at intermediate e, and the phase 
C with the quinoid type structure and the ordinary translational 
period Lf at large e. The primitive unit cells which should be 
used to obtain the translational unit cells of the phases A, B 
and C with the help of rotational (for all phases) and helical 
(only for the phases A and C) symmetry transformations are 
shown in Fig. 1 . 

Since we study the system at constant length and temperature 
T — s> 0, the ground state is determined by the minimum of 
the system internal energy. The Fig. 2(a) shows that at zero 
temperature the phase A is the ground state of the (5,5) nanotube 
for the elongation e < 0.05, the phase B (see Fig. 1(b)) is the 
ground state for 0.05 < e < 0. 13 and the phase C (see Fig. 1(c)) 
is the ground state for e > 0.13. The dependence of the bond 
lengths on the elongation e demonstrates that the structures of 
the phases are different at critical elongation values Sc where 
the internal energies of the phases coincide (see Fig. 2(b)). 
The structure change at £c means that at zero temperature both 
structural phase transitions between the phases A and B and 
the phases B and C are first order transitions with respect 
to the control parameter (elongation). The Young's modules 
of structures are also changed at these phase transitions (see 
Fig. 2(a)). Fig. 2(a) shows the possibility of metastable states 
for the phase A for elongation 0.05 < £ < 0.09, and for the 



phase B for elongation e > 0.13. The calculations show that 
the phase C is unstable for elongation e < 0.13. 

The structure without Stone-Wales (SW) defects is energeti- 
cally favorable for armchair nanotubes at least for elongations 
e < 0.1 [29]. Moreover the high barrier for formation of SW 
defects excludes their spontaneous formation at room temper- 
ature up to the elongation e = 0.15 [29]. Thus we do not 
consider nanotube structures with SW defects. Note also that 
we consider nanotubes at elongations which are well under the 
fracture value (e = 0.24 according to AMI calculations [30] 
and e — 0.3 according to PM3 and density functional theory- 
based calculations [31]). 

Since we are considering a ID system, the term PV in the 
thermodynamic potentials should be replaced with the term 
Fl, where F is the tension force and I is the length of the 
nanotube (a positive force value corresponds to contraction). 
The chemical potential fii for any structural phase i of the 
nanotube takes the form 

fit = Ui- Tsi + FLi, 

where Ui, Si and Li are the internal energy, the entropy and 
the length per one particle of the ID system. It is well known 
that in ID systems at temperature T a co-existence of 
phases with a boundary between them is not possible [32]. The 
phase transitions in ID systems have a crossover character, and 
the crossover region (with respect to the control parameter) in- 
creases with temperature, see Refs. [33,34]. In the crossover 
region the nucleus of another phase is formed as a result of 
thermal fluctuations. From the equality of the chemical poten- 
tials {p,i = fij) of two phases i and j the line of the critical 
values of the tension force Fc and temperature Tc, which cor- 
respond to the middle of the crossover region on the tension 
force-temperature "phase diagram" of the system in equilib- 
rium, can be obtained 

U^-Uj +Tc{Si- Sj) 

r c — . 

Li — Lj 

These structural "phase transitions" in nanotubes can be deter- 
mined by peculiarities of the dependences of elastic constants 
or of the specific heat on temperature or on the tension force. 

The dependence of the internal energy U on the length I of 
the computational cell is interpolated using the Hooke's law 
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where v is the coefficient of elasticity of a CNT with the 
length I. 

The Young's modulus of the (5,5) CNT at the elongation 
e = has the following form 



where R = + -R2)/2 is the average radius of the CNT 
corresponding to the elongation e = 0, w = 0.34 nm is the 
effective thickness of the wall [31,35-37]. As a result the value 
Y K. 1.2 TPa is calculated. The calculated Young's modulus of 
the (5,5) CNT with the Kekule structure is in good agreement 
with the values calculated for the (5,5) CNT with the quinoid 
type structure: 1.1 TPa by PM3 calculations [31], 1.06 TPa by 
Hartree-Fock method [35], 0.95 TPa [36], 0.96 TPa [31] and 
1 .03 TPa [37] by density functional theory-based calculations 
[37]. 

Since the radius R depends on the elongation e, it is not 
convenient to use Eq. (1) for expanded nanotubes. Moreover, 
the coefficient of elasticity v is the quantity which can be mea- 
sured experimentally. Thus, we have calculated the changes of 
the coefficient of elasticity at the structural phase transitions: 
vb — 0.82z^^ and vc = OJii/B, where i^a, vb and vq are the 
coefficients of elasticity for the phases A, B and C, respec- 
tively. 

5. Electronic structure of a deformed nanotube 

In the tight-binding approximation [19], the band structure 
of armchair CNTs with all equal bond lengths has two half- 
filled bands intersecting at the Fermi level E^. They are shown 
by dotted lines in Fig. 3. 

In the case of the quinoid type structure of an armchair 
CNT (Fig. 1(c)) these bands can be expressed as (solid lines in 
Fig. 3) [38]: 

Ei^2{k) = ±ta[l - 2{tb/ta)cos{kLf/2)l 

where ta and ti, are the resonance integrals for bonds a and 
b, respectively, k is the wavenumber of an electron along the 
CNT axis, and is the translational period. In this case the 
armchair CNT remains metallic. 

The calculation of the band structure of an armchair CNT 
with four non-equal bonds (corresponding to the Kekule struc- 
ture of the phase A) is reduced in the tight-binding approxi- 
mation to the diagonalization of the sixth order Hamiltonian 
matrix [10] and the standard zone-folding procedure [19]. The 
standard procedure implies that all atoms are placed in a plane. 
For the primitive unit cell in Fig. 1(a) the Hamiltonian matrix 
has the form: 
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where ta, t^, tc, td are the resonance integrals corresponding 
to the bonds a, b, c, d, respectively; k — kxe^ + kySy is 




2/3 1 

Wavenumber k, n/Li 



Fig. 3. Half-filled bands 1 and 2 intersecting at the Fermi level {E = 0) of 
quinoid type armchair carbon nanotubes (solid lines) for the i.iu, 
tft = to- The dotted lines Ei^2{k) are for equal lengths of bonds (a = b) 
and ta = tt = to. 

the electron wavevector; Bx, By are the unit vectors of the co- 
ordinate system; ri — piqi — ae^, r2 — Piq2 — Psqi — 
-b^e^ - byBy, rs = piq^ = p2qi = -b^e^ + byBy, = 
= P2q2 = de^, rs = p2q3 = -c^Bx - CyBy, rg = 
— —Cx^x + CyBy are the vectors of bonds between atoms. 
For the given (from the PM3 calculations) bond lengths a, b, 
c, d and the translational period Lf we numerically diagonal- 
ized (2) to obtain the electron energy dispersion relation E{k). 
The condition for circumferential quantization kx = 2'iTm/Ch, 
m = —4, . . . , 5, where Ch — 5{2cx + 2d) is the circumference 
length (chiral vector module) of the (5,5) CNT, gives electronic 
energy bands Em{ky) for the nanotube. We assume that the 
dependence of the resonance integral on the C-C bond length 
ac_c is [39]: 

t = [ao/ac-cfto, 

where ag = 0.142 nm, and to = 2.6 eV are parameters for the 
bond length and resonance integral [7,8] of the CNT with all 
equal bond lengths. 

For the undistorted nanotube (all equal bonds) with the ar- 
tificially tripled along nanotube axis translational period, the 
electronic energy band structure in the mapped Brillouin zone 
is shown in Fig. 4(a) for the conduction band near Fermi level. 
In the framework of our approximations (the overlap integral 
for the nearest neighbour atoms is neglected) the valence band 
is the same with the opposite sign. In this case the intersection 
point of half-filled bands is folded to the T point of the Brillouin 
zone. For the distortions shown in Fig. 1(a) the energy (Peierls- 
type) bandgap appears between the valence and the conduction 
bands of the nanotube (see Fig. 4(b)). The value of this gap is 
found to be Eg — 0.24 eV. STM measurements of the electron 
density of states show Eg ~ 0.11 eV for the armchair (7,7) 
CNT at liquid helium temperature [40]. We believe that the 
observed bandgap can be explained by the Peierls transition. 

The diagonalization of matrix (2) for e = 0.029 with the 
PM3 calculated bond lengths of the (5,5) CNT and zone-folding 
gives the electronic energy band structure shown in Fig. 4(c). 
As a result of deformation the lowest conduction band with 
m = is splitted at ky = 0. The diagonalization of (2) for 
fc = gives for the bandgap Eg ~ 2{td ~ tc), i.e., almost 
linear dependence on e according to Fig. 2(b). Note that the 
analogous formula for the bandgap is obtained for the case of 
frani-polyacetylene [19]. 
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Fig. 4. Band structure of the (5,5) carbon nanotube near Fermi level (E = 0): (a) undistorted nanotube (all equal bonds) with artificially tripled period, (b) for 
the case of Kekule ground state distortion (Fig. 1(a)) with a = c = 0.140 nm and b = d = 0.143 nm, and (c) for the phase A at the elongation e = 0.029. 
Only conduction bands are shown, the valence bands are the same with the opposite sign. 



The results of the bandgap Eg calculations for different elon- 
gations e of the phase A are presented in Fig. 5. At the elon- 
gation e = 0.05 the structural phase transition takes place, and 
the bandgap becomes zero for the phase B due to the double 
translational period Lf of this phase (Fig. 1(b)). 
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Fig. 5. Calculated dependence of the (5,5) carbon nanotube energy bandgap 
Eg (at wavevector fc = 0) for the phase A on the elongation e. The vertical 
dotted line indicates the structural phase transition. 



6. Conclusion 

The atomic structure of the (5,5) carbon nanotube is calcu- 
lated for its quasistatic expansion up to the 17% elongation. 
Two first order structural phase transitions with the change of 
the atomic structure symmetry are revealed at the nanotube ex- 
pansion. It is shown that for the elongation greater than 5% 
the structure of the nanotube corresponds at any temperature 
to the metallic phases without the Peierls gap in the electron 
spectrum. The narrow gap semiconductor to metal phase tran- 
sition at nanotube elongation can be used for elaboration of 
nanotube-based stress nanosensors. 
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